Stoner ferromagnetic phase of a graphene in the presence of an in-plane magnetic field 
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We study the effects of an in-plane magnetic field on the ground state properties of both gap- 
less and gapped graphene sheets within Random Phase Approximation. The critical magnetic field 
which leads to a fully spin polarized phase increases by decreasing the carrier density at zero gap 
indicating that no spontaneous magnetic phase transition occurs. However, at large energy gap 
values it decreases by decreasing the density. We find a continuous quantum magnetic phase tran- 
sition (Stoner phase) for Dirac fermions in a doped graphene sheet. Novel in-plane magnetic field 
dependence of the charge and spin susceptibilities are uncovered. 

PACS numbers: 71.10.Ca, 72.10.-d, 73.20.At, 72.25.Dc, 73.50.Fq 



I. INTRODUCTION 

Graphene is an atomically thin two-dimensional (2D) 
electron system composed of carbon atoms on a honey- 
comb lattice. Several experimental groups have recently^ 
introduced techniques which enabled isolation and study 
of systems with one or a small number of graphene lay- 
ers. The interesting physics of graphene systems stems 
from the fact that its envelope-function of low-energy 
Schrodinger equation is equivalent to the massless 2D 
Dirac equation. In the case of graphene the spinor struc- 
ture in the Dirac equation refers to sublattices of its hon- 
eycomb structure and its Brillouin-zone valleys, instead 
of the spin degrees of freedom and the electron-positron 
picture. Graphene, therefore, presents a new type of 
many-body problem in which the noninteracting low en- 
ergy quasiparticle dynamio^ is effectively described by a 
2D massless Dirac Hamiltonian Hq ~ hvpli ■ k, with 
two chiral eigenvalues, ±ek where Ek = /iWF|k|, with 
vp ~ lO^m/sec, is the Fermi velocity of carriers. The 
chirality plays an important role in the novel electronic 
properties of graphene. 

The strength of interaction effects in an ordinary two- 
dimensional electron gas ( 2DEG) increases with de- 
creasing the carrier density. At low densities, the ef- 
fective velocity is suppressed, the charge compressibility 
changes sign from positive to negative, and the spin- 
susceptibility is strongly enhanced. In the Dirac-like 
electrons in graphene, it has been shown^ that the ve- 
locity is enhanced rather than suppressed, and that the 
compressibility, always remains positive, and the spin- 
susceptibility is suppressed. These qualitative differences 
are due to the exchange interactions between electrons 
near the Fermi surface and electrons in the negative en- 
ergy sea. The interband excitations are closely analogous 
to the virtual particle-antiparticle excitations of a truly 
relativistic electron gas. 

Conventional 2DEG has been a fertile source of sur- 
prising new physics for more than four decades. In recent 
years, because of the important and novel physical prop- 
erties found in both theoretical and technological appli- 



cations, there has been a large amount of theoretical and 
experimental studies on the effects of parallel magnetic 
field B in a 2DEG. A great deal of activity was spawned 
in the last decade to understand the apparent metal- 
insulator transition observed in Si-MOSFET and GaAs 
based structuresi^ Although the basic mechanism and the 
existence of a quantum phase transition is still a matter 
of on-going debate, experiments have amassed a wealth 
of data on the transport properties of the 2D electron 
systems in the metallic state. Zhang and Das Sarma^ 
investigated the ground-state properties of the 2DEG in 
the presence of an in-plane magnetic field B using ran- 
dom phase approximation (RPA). They showed that for 
small Wigner-Seitz density parameter = {Trna^)~^^'^ 
in which qb is the Bohr radius in the medium of in- 
terest and in the absence of the magnetic field the sys- 
tem prefers to be in a paramagnetic state. As B in- 
creases the ground-state energy is minimized at a special 
nonzero spin polarization denoted by C*. The degree of 
spin polarization is defined as C = l?^! ~ ^il/('^T + '^i)' 
and J^Kx) is electron density with spin up (down). When 
B increases to a critical value Be in which the system 
is fully spin polarized, there exist two C* values, smaller 
and equal to one, where the total energy is minimized. 
They have shown that the first order phase transition 
from paramagnetic-to-ferromagnetic takes place. Impor- 
tantly, beyond the critical field the energy minimum is 
at ^* = 1 and the system is fully spin polarized. 

Suba§i and Tanatar— , on the other hand, studied the 
same system by using a parameterized expression for the 
correlation energy provided by the quantum Monte Carlo 
(QMC) simulations^ They found that the 2DEG in the 
presence of an in-plane magnetic filed B (B > Be) under- 
goes a first order phase transition to the ferromagnetic 
state (Bloch ferromagnetism) in the density regions as- 
sociated with < Tg < 7 and 20 < < 25, while for 
7 < < 20 their results predicted a continuous phase 
transition (Stoner ferromagnetism). 

In an electron gas system the physical observable quan- 
tities most directly related to the energy are the com- 
pressibility which measures the stiffness of the system 
against changes in the density of electrons and the spin- 
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FIG. 1: (color online). Exchange energy as a function of 
degree of spin polarization, ( for various gap energies. In the 
inset: the correlation energy as a function of ( for various gap 
energies. 
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FIG. 2: (color online). Critical magnetic field as a function 
of inverse square root of density ( in units of 10~^ cm) for 
various gap energies. 



susceptibilities. In bulk electronic systems, the spin- 
susceptibility can usually be extracted successfully from 
total magnetic susceptibility measurements. This is, 
however, likely to be challenging in the case of single-layer 
graphene. In the 2DEG, on the other hand, information 
about the spin-susceptibility can often be extracted from 
weak-field magneto-transport experiments using a tilted 
magnetic field to distinguish spin and orbital response. 

Recently, Hwang and Das Sarma^ have shown that 
the in-plane magnetic field induces graphene magneto- 
resistance which is negative for intrinsic gapless graphene 
while for extrinsic gapless graphene, magneto-resistance 
is a positive value at fields lower than the critical mag- 
netic field and negative above the critical magnetic 
field. The effect of in-plane magnetic field on microwave 
magneto-transporti^ in doped graphene is an open prob- 
lem. 

The purpose of this paper is to study the effects of 
in-plane magnetic field on disorder free doped graphene 
at zero temperature. These effects are very important 
and have some novel and unusual properties in compar- 
ison to the conventional 2DEG. We have revisited the 
problem of the compressibility and spin-susceptibilities 
in the presence of an in-plane magnetic field and find 
that the charge compressibility exhibits a crossover be- 
tween paramagnetic and ferromagnetic phases depending 
on the carrier density and the gap values. Before describ- 
ing the details of the theory and presenting our results, 
we point out that the most novel electronic properties 
of graphene which will be discussed here are based on 
interband interaction and exchange interaction between 
electrons near the Fermi surface in graphene sheets. 

The content of the paper is as follows. In Section II 
we discuss about our theoretical model which contains 
the Zecman energy. Our numerical results are given in 
Section III. Eventually, Section IV contains the summery 



and conclusions. 



II. THEORETICAL MODEL 

We consider a doped graphene sheet with a peculiar 
gap opening due to sublattice symmetry breaking where 
2D massive Dirac fermions at low energy is described by 
noninteracting Hamiltoniani^ Hq = hv-p'a ■ k -I- TOWpcra. 
There are two eigenvalues ±iJk where = y^ej^ 
is the spectrum of particle with A = mwp is the gap en- 
ergy. The results of gapless graphene can be obtained by 
setting A = 0. It should be noted that the in-plane mag- 
netic field couples not only to the spin degrees of freedom 
in a quasi 2DEG which leads to the spin polarization of 
the carriers, but also to the orbital motion of the car- 
riers^'^ 1^^1^^1^'^ due to the finite quantum well thickness. 
However, because of the absence of thickness in graphene, 
the applied B couples only with spin of the carriers and 
leads to the spin polarization of them. Due to the Zee- 
man spin-splitting effect we have a shift in the Fermi wave 
vector for up and down spins k^a = A;f(1 + o'C)^''^ where 
/cp = \/ 4:Tm/gsgv is the unpolarized Fermi wave vector, 
n is the 2D electron density and = 2 and g^ = 2 are 
the spin and valley degeneracy, respectively. The cou- 
pling constant in graphene sheets is density independent 
and given by agr = gsgvB^ /^fivp where e is the average 
dielectric constant of the substrate and air. For a typi- 
cal substrate (e.g. SiC or Si02) the dielectric constant is 
between 1 and 2. 

The total energy per particle in the presence of an in- 
plane magnetic field B as a function of the density ri, the 
spin polarization C, the gap energy A and the coupling 
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constant agr takes the following form 

etot{nX,^,B) = efe„ (n,C, A) +£^(n,C, A) 

+ e,{n,C,A)+SziC,B). (1) 



where 

efei„(n, C, A) 



-^{[ft\|fc|(l + C)+AT/^ 

[fi2i;2fc2(^ _^)_^^2]3/2_2A3} (2) 



is the kinetic energy per particle and 
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is the exchange energy, ^(^^(q, iw, C, A) is the zero tem- 
perature noninteracting polarization function for doped 
graphenoii which is given by 
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where fia- = V h'^v^kpa'^ + A2 and Vq ~ 27re'^/eq is 
the 2D Coulomb interaction. Moreover, the correlation 
energy per particle^ in RPA is given by 



£e(n,C,A) = -£,(n,C,A) + 
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27rn J (27r)2 
In [ ^ — ^)], 



(5) 



and finally the Zeeman energy per particle is ez{C, B) = 
—fisCB where is the Bohr magneton. In the above 
equation we have used the fluctuation-dissipation theo- 
renii^. In order to make the exchange and correlation en- 
ergies finite, we mig ht subtracli^'i^ the vacuum polariza- 
tion energy contributions from the exchange and correla- 
tion energies (S£^(c)(fcF ^ 0) = £j;(c)(A:f) - ex(c)(fcF = 0). 
Due to the conserved number of states in the Brillouin 
zone, we do need an ultraviolet momentum cut-off 
which is approximated by 7rfc2 = (27r)2/y^g where is 
the area of the unit cell in the honeycomb lattice. The 
dimensionless parameter A is defined as kc/kp. 

The total energy per particle for a gapless graphene in 
the noninteracting electron scheme is given by 



.0 



t(.n,(:,B) 



[(l+C)'/' + (l-C)'/']-A'i3CS (6) 
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where £f = hv-pk-p is the Fermi energy of the gapless 
graphene. The minimum of the noninteracting energy as 



a function of spin polarization occurs at Cq = 2^BB{ep — 
^^i?2)i/2/£2 ^ Setting = 1 allows us to determine the 
critical magnetic field Bco{n) corresponding to the fully 
spin polarize the system. The critical magnetic field for 
the noninteracting system is Bco = £f/ V^/^ b ■ 

To calculate C*{B) for the interacting case, the total 
energy in Eq. ([1]) needs to be minimized with respect to 
C and then the critical magnetic field Be can be found 
for the fully spin polarized case. At a finite magnetic 
field the energy minimum occurs at nonzero polarization 
< < 1- Beyond the critical field the energy minimum 
is at C* = 1 and the system is fully polarized. In general, 
the critical magnetic field takes the form 



B^ 

Bco 



2eF 



(2£|, + A2)1/2_ A]+2 



IC=i 



(7) 



III. NUMERICAL RESULTS 



We now turn to the presentation of our numerical re- 
sults. We consider agr = 1 which is a appropriate value 
for graphene placed on the SiC substrates and we choose 
the gap values between and 100 meV which is observed 
in typical experiments. 

In Fig. 1 we plot the exchange ( in the inset: correla- 
tion ) energy of graphene as a function of ^ for a range 
of A values. The exchange energy is positive because our 
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FIG. 3; (color online). Total energy as a function of spin 
polarization for various magnetic fields for (a): A = and 
(b): A = lOOmeV at A = 100. 



regularization procedure implicitly selects the chemical 
potential of undoped graphene as the zero of energy. It 
would be noted that d5ex/dC,\c,=i > at small A and it 
changes the sign at large gap values. The slope of ex- 
change and correlation energies with respect to C around 
C = 1 have opposite signs. Note that both fex and |&c| 
have the same density dependence and they both increase 
with decreasing density^ These arguments will be use- 
ful in describing the the critical magnetic field given by 
Eq.El 

In Fig. 2, we plot the calculated critical field Be which 
polarizes the quasiparticles for the interacting case in the 
units of the critical field for noninteracting Dirac mass- 
less fermions Bcq, as a function of inverse square root of 
the density. The critical field increases by decreasing the 
density of carriers for massless case due to the impact of 
exchange energy ( see Eq. [71) . This particular feature is 
in contrast to the 2DEG in which the reduction of car- 
rier density leads to the decline of the critical fieldi^iii 
This distinguished behavior is a direct result of the chi- 
rality in the massless Dirac fermions and interband in- 
teraction features. The critical fields, on the other hand, 
have a non monotonic behavior at small gap values. It 



increases by decreasing the density till it reaches to the 
maximum value and then decreases at very low density 
values because of the competition between the exchange 
and correlation energy contributions, as shown in Fig. 1. 



We have found that Br vanishes at about 
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for A = 100 meV. This is a similar behavior to the con- 
ventional 2DEG where system goes to the fully polarized 
state spontaneously at a special density, Vg ~ 25.5 calcu- 
lated within Monte Carlo simulations^ or Ts k, 5.5 base 
on RPA calculation.— 

Moreover, at low density region associated with large 
A, we found no indication for a spontaneous magnetic 
phase transition for small A even at large coupling con- 
stants. These results are in contrast to the results re- 
ported in Ref. [l^ where the exchange term, the Hartree- 
Fock theory was only used. In that work, the authors 
found that exchange interactions between Dirac fermions 
can stabilize a ferromagnetic phase at low doping when 
the coupling is sufficiently large. We have not found any 
evidence for this instability using RPA calculations. The 
RPA is a minimal dielectric scheme that allows quanti- 
tative predictions beyond the Hartree-Fock theory. In 
the present case of a two-dimensional electron gas on a 
graphene sheet, the Hartree-Fock exchange contribution 
to the ground-state energy is positive. In our work we 
clearly show that the RPA correlation energy is negative. 

Furthermore, it is showoi^ that the kinetic energy en- 
hancement of the spin-polarization phase nearly cancels 
the exchange enhancement and the correlation energy 
plays a dominant residual role. Therefore, the inclusion 
of the correlation energy suppresses the spin-polarized 
phase found in the exchange only calculation in gapless 
graphene. In examining the tendency of the system to de- 
velop magnetic order in the presence of electron-electron 
interactions it is thus crucially important to include both 
exchange and correlation contributions. 

In Fig. 3 we plot the calculated total ground state en- 
ergy, fctot(fcF) = Etotik-p) - Stotikp = 0) in units of the 
Fermi energy for massless Dirac fermions e-p as a func- 
tion of the spin polarization parameter ^. The results 
are shown for various magnetic fields at (a) A = and 
(b)A = lOOmeV. In both cases the minimum energy oc- 
curs at paramagnetic state, namely C* = in the ab- 
sence of magnetic field but as B increases the minimum 
energy shifts to non-zero spin polarization and C* in- 
creases continuously to the ferromagnetic phase (C* = 1) 
at B = Be- For B > Be the system remains in the ferro- 
magnetic phase. This indicates a continuous-phase tran- 
sition (Stoner type) from para- to- ferromagnetic phase 
in the presence of magnetic field for each density value 
whereas a first-order phase transition for whole of the 
density range is predicted for 2DEGJ^ 

In Fig. 4 we plot the magnetization C*(^) as a func- 
tion of the applied magnetic field BiLiLL Clearly there is 
no longer jump in the magnetization ai B = Be due to a 
continuous-phase transition in graphene. In the conven- 
tional 2DEG, the transition to the ferromagnetic state 
near the critical magnetic field value happens with a dis- 
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FIG. 4: (color online). Spin polarization as a function of the FIG. 5: (color online). Compressibility of gapless graphene as 

magnetic field for several energy gap values at A = 100. The a function of inverse square root of density (in units of 10~® 

inset: The spin susceptibility as a function of the magnetic cm) for both fully spin polarized and unpolarized states, 
field for A = 0. 



Crete jump in the polarization indicating a first order 
transition to the fully polarized state. The magnetiza- 
tion (^* is a semi-linear function of the magnetic field 
versus the large gap values. 

A quantity of interest which can be accessed experi- 
mentally is the non-linear spin susceptibility of the sys- 
tem defined as x/xoiB = 0) = j^dC* /dB where xo is 
the Pauli susceptibility. The spin susceptibility decreases 
nonlinearly by increasing the magnetic field at small A 
values, showing that the polarizability of the system de- 
creases. This feature should be verified by magneto- 
resistance measurements through the polarization field 
Be 

Another important thermodynamic quantity is the 
compressibility, k which yields interesting features when 
graphene is subjected to an in-plane magnetic field. The 
exchange energy is positive while the correlation energy 
is negative. This has important implications on the ther- 
modynamic properties. The compressibility can be cal- 
culated from its definition, = n'^d'^{nSetot)/dn'^. In 
Fig. 5 we have shown the inverse of compressibility of 
gapless graphene as a function of the inverse square root 
of density for unpolarized and fully polarized states, kq/k, 
increases with decreasing density at small gap energy. 
This behavior is in contrast to the conventional 2DEG. 
The compressibility of noninteracting gapless graphene 
is kq = 2/(n£F)- The exchange energy tends to reduce 
the compressibility while correlations tend to enhance it. 
We found that at the given Bc{A = 5) which is associ- 
ated with a special density, kq/k starts from C = 1 ^^'^ 
slowly tends to the paramagnetic results of the gapless 
one. This special behavior is a consequence of the fact 
that Bc{A < 5) is smaller than i?c(A = 5) as shown in 
Fig. 2. However at gapped graphene, say A = 100, we 
chose a value of the critical magnetic field B^A = 400) 
and observed that kq/k switches to its fully polarized 



system value with a kink-like behavior as shown in the 
Fig. 6. This feature is a consequence of the fact that 
Sc(A < 400) is larger than B^A = 400). This suggests 
that in the compressibility measurements the effect of the 
polarizing magnetic filed could be discerned. The phys- 
ical reason for having two different behaviors at small 
and large A is that the critical magnetic field behaves 
in different ways at small and large energy gap values. 
Note that at very large gap energy, kq/k decreases by 
increasing The non-monotonic behaviors of kq/k 

with respect to A is due to the comparison between the 
exchange energy and the correlation energy as a function 
of gap values. 



IV. CONCLUSION 

In summary, we study the effects of in-plane magnetic 
field on the ground state properties of both gapless and 
gapped graphene where the conduction band is partially 
occupied. The present work is demonstrative the increas- 
ing behavior of the critical field in which the system be- 
comes fully spin polarized by decreasing the density for 
gapless graphene. Accordingly, it means that there is no 
longer a spontaneous para-to- ferromagnetic phase tran- 
sition for gapless graphene at zero-magnetic field. The 
critical magnetic field decreases by decreasing the den- 
sity at large gap values. Quite interestingly, we find a 
continuous quantum magnetic phase transition for the 
whole range of the density at zero temperature. The 
novel in-plane magnetic field dependence of charge and 
spin susceptibilities are obtained. The inverse compress- 
ibility as a function of inverse density exhibits a crossover 
from the fully polarized state to the paramagnetic case 
for gapless graphene, which should be identifiable exper- 
imentally. 
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FIG. 6; (color online). Compressibility of gapped graphene 
with A = lOOmeV as a function of inverse square root of 
density (in units of 10"^ cm) for both fully spin polarized 
and unpolarized states. 



It is worth noticing that the aforementioned results are 
in contrast to those results calculated for a conventional 
2DEG due to the effect of interband interaction. We 
suggest these results should be verified by experiments. 

It is recently shown that ripples in graphene induced 
a gauge field. It is convenient to emphasize that the 
study of the effects of parallel magnetic field on physical 
quantities in graphene sheets at the presence of such an 
induced gauge filed is an interesting problem which might 
be taken into account. 
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